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Abstract 


The  review  of  the  literature  on  this  subject  has  been  accomplished. 
The  equations  and  computer  programs  for  the  conversion  of  energy  have 
been  written  and  tested  on  the  computer.  This  research  has  shown  that 
the  hemispheric  meteorological  data  from  the  AFGWC  and  from  ETAC  have 
serious  gaps  that  hinder  the  progress  of  this  task.  For  future  work  it 
is  suggested  to  introduce  an  interpolation  technique  which  will  formally 
fill  the  data  gaps,  or,  alternatively,  to  acquire  new  data  sets  without 

such  gaps.  Then  the  calculation  of  energy  transfers  can  be  accomplished. 
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The  equations  describing  atmospheric  motion  form  a  nonlinear  system. 
The  motions  range  from  the  planetary  scale  to  the  microscale.  The  kinetic 
energy  of  large-scale  atmospheric  motion  is  transferred  between  the  scales 
of  motion  by  the  nonlinear  interaction  of  atmospheric  waves.  The  purpose 
of  this  research  is  to  study  two  formulations  of  the  nonlinear  transfer  of 
kinetic  energy  between  the  larger  scales  of  motion  using  a  data  set  which 
only  recently  has  become  available  and  which  has  not  been  used  for  this 


purpose 


2 


Many  studies  of  kinetic  energy  transfer  have  concentrated  on  the 
total  eddy  kinetic  energy  by  utilising  the  governing  equations  in  a 
physical  space  formulation.  A  great  deal  also  can  be  learned  by  look¬ 
ing  at  the  growth  and  decay  of  the  Individual  scales  of  eddies.  It  is 
thus  necessary  to  convert  the  governing  equations  to  spectral  fore  or 
to  the  domain  of  scale  or  wave  number.  Two  methods,  which  will  be 
described  in  Section  3,  have  been  proposed  to  study  the  nonlinear  trans¬ 
fer  of  kinetic  energy  in  the  spectral  domain. 

Spectral  models  increasingly  are  being  tested  for  use  In  numerical, 
weather  prediction  (Machenhauer  and  Daley*  1972;  Bourke,  1974;  Bourke 
et  al. ,  1977)  and  extended  range  prediction  (Baede  and  Hansen*  1977) . 
Recent  studies  with  general  circulation  models  have  examined  the  model¬ 
generated  kinetic  energy  using  spectral  analyses  (Stone  et  al. *  1977; 
Baker  et  al.,  1977*  1978).  Observational  spectral  studies  of  atmospher¬ 
ic  kinetic  energy  increasingly  are  needed  on  at  least  a  hemispheric 
scale  for  monitoring  and  evaluating  the  efficiency  of  numerical  weather 
prediction  models  and  general  circulation  models.  Despite  the  work  of 
many  researchers*  fully  hemispheric  spectral  statistics*  and  to  a  lesser 
extent  their  analysis  and  interpretation,  are  still  only  poorly  known. 

Lack  of  a  single  hemispheric  analysis  scheme  has  hampered  previous 
studies  of  the  nonlinear  transfer  of  kinetic  energy.  Past  researchers 
used  only  a  middle-latitude  analysis,  a  tropical  analysis,  or  meshed  the 
two  analyses  to  obtain  fully  hemispheric  data  coverage.  Separate  middle- 
latitude  and  tropical  analyses  overlap  in  areal  extent  but  do  not  agree. 
Meshing  of  the  two  analyses  could  Introduce  extraneous  waves  which  would 
contaminate  the  statistics  on  the  nonlinear  transfer  of  kinetic  energy. 
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The  National  Meteorological  Center,  in  1974,  and  the  Air  Force  Glo¬ 
bal  Weather  Center,  in  1976,  began  using  the  Flattery  analyses  schema 
(Flattery,  1970;  Bergman  et  el.,  1974).  This  analysis  scheas  is  very 
different  from  the  previous  analysis  scheaes  and  gives  one  analysis  for 
the  entire  hemisphere.  This  research  will  use  data  from  the  Flattery 
analysis  scheme  for  computation  of  the  nonlinear  transfer  of  kinetic  ■ 
energy. 

2.  Present  status  of  the  question 

Two  methods  have  been  proposed  for  the  study  of  energetics  in  the 
domain  of  wave  number.  Saltrman  (1957)  proposed  a  one-dimensional  wave 
number  formulation  and  derived  energy  equations  which  govern  the  beha¬ 
vior  of  the  individual  scales.  Dependent  variables  along  a  latitude 
circle  are  expressed  in  terms  of  a  Fourier  series.  The  wave  nusber 
represents  the  number  of  waves  around  the  latitude  circle. 

Baer  (1974)  proposed  the  degree  of  the  associated  Legendre  poly¬ 
nomial  as  a  two-dimensional  spectral  index  where  the  dependent  varia¬ 
bles  are  expressed  in  terms  of  spherical  harmonics.  This  index  incor¬ 
porates  two  one-dimensional  indices.  Burrows  (1975)  used  this  approach 
to  derive  kinetic  energy  equations  for  the  individual  scales. 

The  equations  of  horizontal  motion  and  the  continuity  equation  in 
spherical  coordinates  sre 

$  *<u5q> &+  £  J&- fv  + 


la  these  equations  tha  Independent  variables  (X,$,p,t)  represent  lon¬ 
gitude,  latitude,  pressure,  and  tine,  respectively.  The  dependent 
variables  are  u,  v,  and  a,  the  west  and  south  wind  eonponents,  and  the 
geopotential  height  of  the  lsobarlc  surface,  respectively;  a  is  the 
radius  of  tha  earth;  f  -  2Qsin+  is  the  coriolis  parameter ,  where  Q  is 
the  earth's  angular  velocity;  g  la  tha  gravitational  acceleration;  and 
X  and  T  are  the  eastward  and  northward  eonponents,  respectively,  of  the 
frictional  force  per  unit  nase.  These  equations  are  used  to  derive  the 
kinetic  energy  equations  in  the  one-  and  two-dimensional  formulations 
described  In  the  following  paragraphs.  The  above  equations  do  not  In¬ 
clude  the  vertical  advectlon  terns.  The  reasons  for  this  are  explained 
In  Section  4. 

A.  One-dimensional  wave  number  formulation 

The  basis  for  Salt  man's  formulation  Is  that  any  of  the  meteorolo¬ 
gical  dependent  variables  specified  along  a  latitude  circle  may  be 
written  in  terms  of  a  Fourier  representation 

fW-1 
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The  complex  coefficients,  F(t),  are  given  by 

and  are  the  representation  of  f (X)  In  the  wave  number  domain. 

The  spectral  kinetic  energy  equation  originally  given  by  Saltsman  ’ 
(1957),  has  been  modified  (e.g.,  Tang,  1967)  to  apply  to  a  limited 
latitudinal  domain  and  Is 

W<P*P= 
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where  *^(£)*[F(J.)G(-l)+F(-A)G(£)]  represents  the  correlation  of  two 
functions,  and  <jifg(m,l)«*[F(i-m)G(-t)+F(-i-m)G(l)]  represents  part  of 
the  triple  correlation  product.  K(l)  is  the  kinetic  energy  in  wave 
number  1;  a  is  again  the  radius  of  the  earth;  u  and  v  are  the  mean 
zonal  and  meridional  wind  speeds.  Independent  variables  as  sub-scripts 
show  differentialtlon  with  respect  to  that  variable. 

The  transform  pairs  are  f(l):  u  v  z 

F(t)s  U  V  A. 

The  first  integral  on  the  right-hand-side  of  the  equation  repre¬ 
sents  the  conversion  from  the  zonal  kinetic  energy  into  the  kinetic 
energy  of  the  1th  component  through  the  work  of  Reynolds  stresses  against 
the  zonal  velocity  gradients.  This  is  called  a  wave-zonal  interaction. 

The  second  Integral  represents  the  transfer  of  kinetic  energy  of 
other  wave  components  into  that  of  the  1th  component  through  the  non¬ 
linear  wave  interaction.  Summation  of  this  term  over  all  wave  numbers 
is  zero  since  the  nonlinear  interaction  does  not  create  or  destroy  energy, 
but  only  redistributes  it  between  the  scales  of  motion.  This  is  called  a 
wave-wave  interaction. 


The  third  and  fourth  terms  are  fluxes.  The  third  term  la  the  flux 
of  eddy  kinetic  energy  arising  from  the  nonlinear  Interaction  among  the 
waves.  The  fourth  term  represents  the  work  done  by  the  pressure  force 
at  the  boundaries . 

The  symbol  ♦  .  (i)  also  can  be  written  as  l  ♦f_(l)“fTfP\ 
r*  n«l  r® 

Accordingly,  the  expression  for  the  spectrum  of  the  meridional  transport 
of  any  quantity  f  is  given  by 

Several  authors  have  studied  the  effects  of  specific  terms  In  the 
above  formulation  on  the  maintenance  of  the  kinetic  energy  of  the  indi¬ 
vidual  scales.  Saltzman  (1970)  reviewed  previous  studies.  There  are 
two  primary  studies  dealing  with  the  wave-wave  and  wave-zonal  interac¬ 
tions.  Both  of  these  studies  used  only  the  horizontal  wind  components 
in  middle-latitude  regions.  The  general  results  from  these  studies  will 
be  discussed  below. 

Saltzman  and  Teweles  (1964)  used  nine  years  of  500  mb  data  to  derive 
statistics  on  the  rate  of  transfer  of  kinetic  energy  in  the  wave-wave  and 
wave-zonal  interactions.  They  used  geostrophlc  winds  and  Included  up  to 
wave  number  t-15.  In  the  mean,  kinetic  energy  is  transferred  from  the 
eddies  to  the  mean  flow  with  maximum  contribution  from  wave  numbers  2  and 
7.  The  minimum  contribution  is  from  wave  number  4. 

Wave  nuofoers  1,  3,  4,  and  11-15  gain  kinetic  energy  from  the  wave- 
wave  Interaction  with  wave  numbers  2  and  5-10.  Saltzman  and  Teweles 
(1964)  first  suggested  that  the  loss  by  wave  number  2  is  replenished  by 
a  large  forced  conversion  of  available  potential  energy  on  this  scale 
associated  with  the  ocean-continent  structure.  The  loss  by  wave  numbers 
5-10  is  replenished  by  a  conversion  from  available  potential  energy  result 
ing  from  baroclinic  instability.  This  wave  number  band  is  the  one  for 
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whlch  baroclinic  instability  is  the  largest  and  corresponds  to  the  cyclone- 
scale  waves  to  both  larger  and  smaller  waves  is  in  general  agreement  with 
the  theoretical  results  of  Fjj/rtoft  (1953). 

Yang  (1967)  used  one  year  of  data  at  eight  pressure  levels  from  lOOOmb 
through  100  mb  to  study  the  vertical  variations  of  the  interactions.  He 
showed  that  the  conversion  to  zonal  kinetic  energy  reaches  a  maximum  at  the 
jet  stream  level.  The  contribution  of  the  long  waves  increases  with  height 
to  the  jet  stream  level  whereas  the  contribution  from  the  short  waves  reach¬ 
es  a  maximum  slightly  below  jet  stream  level.  For  the  wave-wave  inter¬ 
action  at  low  levels,  the  medium  and  short  waves  provide  kinetic  energy  to 
the  long  waves.  However,  throughout  most  of  the  troposphere  the  medium 
waves  provide  kinetic  energy  to  both  long  and  short  waves. 

The  previously  discussed  results  apply  to  middle-latitudes  and  are 
much  different  from  the  interactions  in  the  tropics.  Kanamitsu  et  al. , 

(1972)  concluded  from  200  mb  data  that  the  ultralong  waves  1*1,2  are  im¬ 
portant  in  the  kinetic  energy  budget  for  the  tropics.  A  thermally  direct 
east-west  vertical  circulation  (Krishnamurti,  1971)  due  to  land-ocean 
heating  contrasts  evidently  generates  kinetic  energy  in  wave  numbers  1  and  2. 

For  wave-zonal  Interactions  in  the  tropics  1*1  provides  the  greatest 
amount  of  energy.  All  waves  1  >.  **  draw  energy  from  the  zonal  flow  with  a 
peak  at  about  wave  number  8.  For  wave-wave  interactions  1*1,2  lose  energy 
to  all  other  waves. 

Saltzman  and  Flelsher  (1960)  noted  that  the  same  wave  number  may 
actually  represent  vastly  different  physical  scales  at  widely  separeated 
latitudes.  The  two-dimensional  spectral  index  formulation,  to  be  dis¬ 
cussed  next,  attempts  to  overcome  this  ambiguity  and  to  more  adequately 
represent  the  two-dimensional  nature  of  the  large-scale  atmospheric  flow. 


B.  Two-dimensional  spectral  index  formulation 

The  surface  spherical  harmonics,  Y£q  (X,y)«Pln  (y)exp(i£X) ,  are  also 
appropriate  functions  for  a  spectral  expansion  of  meteorological  varia¬ 
bles.  The  series  expansion  of  a  scalar  variable  can  be  written  as 

AU,M)=  t  f 

-U-L  n=Ui 

The  coefficients  in  this  expansion  are  given  by 

~  HT  J.  f  A(XjM)  (XjM.)  Jlx 

-J  o 

The  associated  Legendre  polynomials  are  given  by 


In  the  above  equations  y»sin$**sin (latitude) ;  L  and  N  are  the  trunca¬ 
tion  planetary  wave  number  and  planetary  index,  respectively;  Y*  (X,y)“P.  (p) 

In 

exp  (iiX)  is  the  complex  conjugate  of  Y^CX.y) ;  and  P_fcn(y)=P&n(y)  (Platzman, 
1962).  For  the  spherical  harmonics,  n-i  gives  the  number  of  zeroes  between 
the  poles.  Baer  (1972)  proposed  the  use  of  the  degree,  n,  of  the  associated 
Legendre  polynomial  as  a  two-dimensional  scale  index  or  planetary  index. 

The  full  set  of  expansion  coefficients  can  be  found  only  when  the 
meteorological  variables  are  known  over  both  hemispheres.  Since  meteoro¬ 
logical  variables  for  only  half  the  atmosphere  will  be  used,  it  is  necessary 
to  make  the  further  assumption  that  the  fields  are  either  symmetric  (even 
parity  expansion)  or  antisymmetric  (odd  parity  expansion)  about  the  equator. 
Symmetry  is  defined  by  n-£  even  and  antisymmetry  by  n-l  odd. 

Burrows  (1975,  1976)  divided  the  wind  components  into  rotational  (R) 
and  divergent  (D)  parts  given  by  u*uR(even)+Uj>(odd) ,  v*vR(odd)+Vp(even) 
where  the  even/odd  refer  to  the  symmetry/antisymmetry.  The  procedure  for 
separation  of  the  wind  components  is  to  find  the  rotational  part  of  the 
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wind  and  then  subtract  the  wind  in  grid  points  to  find  the  divergent  com¬ 
ponents  Ujj^u-u^,  Vjj*v-v^.  The  divergent  components  are  then  expanded  in 
spectral  series  from  the  grid  point  values. 

The  procedure  for  finding  the  rotational  part  of  the  wind  field  de¬ 
pends  on  a  decomposition  of  the  wind  field  into  a  sum  of  irrotational  and 
nondivergent  parts  (e.g.,  Holton,  1972).  The  vorticity  is  found  at  each 

grid  point  and  expanded  in  a  series  to  find  the  coefficients  .  The 

-a2 

coefficients  of  the  stream  field  1 1/„  are  found  from  =  —s  ■  where 

r£n  in  n(n+l)  in 

use  is  made  of  the  identity  V^Y.  =*  n-^r^  Y  .  The  stream  field  is  then 

xn  z  Jen 
a 

found  at  each  grid  point  from  the  spectral  series.  Rotational  wind  com¬ 
ponents  are  found  using  finite  differences  and  the  relations 


-/" (1-u2)  3<l» 

V  —a  -  ;  V" 


*  *R  a/" (l-pz) 

Finally,  series  are  found  for  u^  and  vR  from  the  grid  point  values. 

The  spectral  kinetic  energy  equations  were  given  by  Burrows  (1975). 
The  symbolic  formula  for  the  kinetic  energy  of  one  spectral  wind  compo¬ 
nent  is 

=  (i-  e*  ^ 

2Ai 

Ol  “  *  ci  Tju 

MAI  MAz 

+  ^  ^  zxlmvz  -4=^4* 

S’PMi  SW>-  Co*  P.<3. 

♦  F«J  e_t<a  P4n  izt  . 

FAlc  i+n  even 


These  equations  are  symbolic  because  it  is  understood  that  series  are  to 
be  inserted  for  wind  components  and  height.  For  example,  the  MAI  tern 
in  the  above  equation  would  be  written 


J\+n,  even 


x  odd 


Ai  4-n3  odd 


There  are  similar  equations  for  )J  »atid 

)J-  1°  the  above  equations,  the  symbols  below  the  terms  stand  for 

zonal  and  meridional  advection,  sphericity  terms,  coriolis,  pressure  gra¬ 
dient,  and  "friction  remainder",  respectively.  5on={l*  n=0  ^  is  t*le 

Kronecker  delta.  The  nonlinear  interaction  among  the  scales  of  motion  is 
given  by  the  sum  of  advection  and  spfericity  terms. 


The  pressure  gradient  term,  in  the  above  equations,  does  not  generate 
kinetic  energy.  It  does  act  to  balance  wind  and  mass  fields  which  are 
thrown  out  of  balance  by  turbulence  and  generation  of  kinetic  energy  from 
available  potential  energy.  The  generation  of  kinetic  energy  from  avail¬ 
able  potential  energy  is  not  present  in  these  equations  because  of  the 
assumed  parities  for  the  wind  and  height  fields.  The  pressure  gradient 


force  contributes  only  to  the  rotational  wind  components  which  are  non- 
divergent  and  thus  cannot  show  the  generation  of  kinetic  energy  from 
available  potential  energy. 


The  coriolis  force  does  no  net  work  but  affects  the  total  kinetic 
energy  in  each  n-scale.  It  should  be  noted  that  the  coriolis  term  does 
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no  t  enter  into  the  equations  based  on  grid  point  values.  It  does  enter 
into  the  i-domain  (one-dimensional  wave  number)  equations  for  the  indi¬ 
vidual  components,  but  when  the  equations  are  added  to  give  the  total 
kinetic  energy  in  each  i-scale,  the  corlolls  contribution  is  zero. 

Burrows  (1975,  1976)  has  made  the  only  study  of  the  nonlinear  trans¬ 
fer  of  kinetic  energy  using  the  two-dimensional  formulation.  He  used  15 
days  of  August,  1970  data  at  eight  pressure  levels.  He  also  compared  his 
n-domain  (two-dimensional  index)  results  to  the  averages  of  Saltzman  and 
Fleisher  (1962)  for  middle-latitudes  and  to  the  results  of  Kanamitsu  et  al. , 
(1972)  for  tropical  latitudes.  It  should  be  noted  that  sho&term  averages 
of  the  nonlinear  transfers  may  not  be  representative  of  longer  term  aver¬ 
ages  . 


3.  Prepared  computer  programs 

The  primary  contribution  of  this  research  is  an  application  of  Burrow's 
two-dimensional  formulation  for  the  nonlinear  transfer  of  kinetic  energy  to 
a  longer  period  in  a  different  year  and  involving  objectively-analyzed  data 
on  a  fully  hemispheric  grid.  Saltzman's  procedure  is  applied  to  the  same 
data. 

Data  for  August  1977  have  been  obtained  for  the  whole  Northern 
Hemisphere  from  AFGWC  on  magnetic  tapes  and,  independently,  from  ETAC.  The 
data,  obtained  by  the  Flattery  scheme,  are  for  the  points  on  a  numerical 
weather  prediction  grid  extending  to  the  equator  at  eight  pressure  levels 
from  1000  mb  through  100  mb.  For  the  purpose  of  using  the  present  programs, 
the  data  are  interpolated  to  a  latitude/longitude  grid  by  interpolation, 
using  the  interpolating  polynomial  through  12  nearest  points. 
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After  Interpolation  to  a  latitude/longitude  grid,  the  wind  and  height 
for  each  pressure  level  are  converted  to  the  spectral  fora  by  the  methods 
described  in  Section  2.  The  spectral  coefficients  are  the  basic  variables 
for  the  remainder  of  the  calculations • 

For  an  Initial  analysis  of  the  nonlinear  transfer  of  kinetic  energy 
based  on  data  from  the  new  analysis  scheme  it  is  advantageous  to  retain 
only  those  terms  of  major  importance.  The  net  effect  of  the  nonlinear 
terms  involving  vertical  velocity  is  one  order  of  magnitude  less  than  the 
major  horizontal  terms.  In  addition,  vertical  velocity  is  difficult  to 
evaluate  even  with  a  large  expenditure  of  resources.  The  terms  involving 
vertical  velocity  are  therefore  excluded  from  the  kinetic  energy  equations. 

Each  term  in  the  kinetic  energy  equations  is  evaluated  for  each  pres¬ 
sure  level.  The  time  derivatives,  from  the  left-hand-side  of  the  equations, 
is  evaluated  by  centered  finite  time  differences.  The  terms  from  the  right- 
hand  side  of  the  equations  are  evaluated  from  values  of  height  and  wind  at 
a  single  synoptic  time,  with  the  exception  of  the  "friction  remainder"  term 
which  is  used  to  balance  the  equation.  The  "friction  remainder"  contains 
not  only  friction  but  also  the  contribution  from  unevaluated  terms,  errors 
in  the  evaluations,  and  errors  arising  from  equating  time-centered  and 
instantaneous  values.  The  unevaluated  terms  are  the  generation  of  kinetic 
energy  by  the  conversion  of  available  potential  energy,  the  interaction 
with  waves  smaller  than  the  resolved  scale,  and  in  the  two-dimensional 
formulation  the  fluxes  of  kinetic  energy  due  to  cross-equatorial  flow. 
Unresolved  scales  also  influence  the  kinetic  energy  in  the  resolved  scales 
due  to  aliasing. 

The  nonlinear  terms  in  the  two-dimensional  formulation  are  evaluated 
by  the  transform  method  described  by  Machenhauer  and  Rasmussen  (1972) 
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whlch  eliminates  the  requirement  for  large  numbers  of  interaction  coeffi¬ 
cients  (Platzman,  1960).  Very  briefly,  the  transform  method  involves  a 
conversion  from  the  spectral  domain  to  the  physical  domain,  where  the  mul¬ 
tiplication  is  performed,  followed  by  a  conversion  back  to  the  spectral 
domain.  A  Fast  Fourier  Transfrom  program  written  by  Singleton  (1969)  is 
modified  and  used  for  integration  over  longitude.  Gauss-Legendre  quad¬ 
rature  is  used  for  integration  over  latitude. 

A  listing  of  all  these  programs  is  available  from  Major  Frank  Bower, 
AFGWC. 

The  results  from  the  computation  are  not  available  since  the  data  on 
the  tapes  from  the  National  Weather  Service  and  from  ETAC  have  serious 
gaps.  The  transfrom  methods,  on  the  other  hand,  give  reliable  results 
only  with  complete  data  sets,  without  gaps.  For  this  reason  there  will 
be  a  delay  of  a  year  or  so  until  we  find  a  suitable  interpolation  method 
or  some  other,  more  complete,  data  sets. 
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